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We conduct a careful analysis of the data provided bv lKrogstad fe Davidson ( 2011 ) and 
show that their data do not support their conclusions. According to their published 
data, their decaying approximately homogeneous isotropic turbulent flows are, invari- 
ably, clearly different from Saffman turbulence; and very clearly marked differences exist 
between the far downstream turbulence behaviors generated by their conventional grid 
and by their multiscalc cross grids. 



1. Introduction 



A few years ago, lLavoie. Dienidi fe Antonial (|2007|) investigated potential effects of 
initial conditions on the decay of approximately homogeneous isotropic turbulence. Initial 
conditions refer to the way the turbulence is generated. In the wind tunnel experiments 
of these authors, the turbulence was passively generated by square-mesh biplane grids 
placed at the test section entry. A particular aspect of the potential dependence on initial 
conditions is whether the power-law decay of the far-downstream turbulence depends on 
them. Quantitatively, the question is whether the decay exponent n in 



(x - x y 



(i.i) 



(where u 2 stands for two thirds of the turbulent kinetic energy and x is the streamwise 
distance alon g the tunnel, xp being a virtual origin) differs for different initial conditions 
as claimed by George! ( 19921) . 

Lavoie et alT l 20071 ) tried four different conventional passive grids (with square or with 



round bars with/ without a small helical wire) and two different test sections (one with 
and one without a secondary contraction to improve isotropy) . They did not find any sig- 
nificant effect of initial conditions on the decay exponent n other than that of anisotropy 
w hich does, itself, depend on i nitial conditions and persists far downstream. 



Krogstad fe Davidson! (|2011l ) carried out a similar wind tunnel study but with two mul- 
tiscalc grids and one conventional grid. Their grids were all monoplanar and their two 
m ultiscale grids were chosen from one of the three families of multiscalc grids introduced 



bv lHurst fc Vassilicosl ( 2007 ). specifically the family of fractal cross grids. These grids are 



very different from the low-blockage space-filling fractal square grids which have been the 



generated turbulence (ISeoud & Vassilicos 2007: 


Nagata. Suzuki. Sakai. Havase & Kubo 


2008all&: Stresing. Peinke. Seor 


Ld & Vassilicosl201Ci 


Mazellier & Vassilicosl2010l:ISuzuki. Nagata. Sakai & Rvota 


20101 Laizet & Vassilicosl2011; 


Valente & Vassilicos 


|2011 ). The reason whv multiscale /fractal 



cross grids have mostly been negl ected (except in studies where they were used to 
enhan c e the Reynolds number, see Kinzel. Wolf. Holzner. Liithi. Tropea fc Kinzelbachl 



( 2010l ); lGeipel. Henry Goh fc Lindstedtl (|2010l )) is that lHurst fc Vassilicosl (|2007t ) did not 
make any strong or unexpected claim about the dependence of u 2 on x — xq in decaying 
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turbulence generated by them. Their conclusion on these grids was just a double nega- 
tive: "the turbulence decay observed is not in disagreement with power-law fits and the 
principle of large eddies" . 

Sketches of the multiscale cross grids used by Krogstad fe Davidson! ( 2011 ) can be seen 
in their figure 1 and are described in their section 2 where they are labeled msgl and 
msg2. We do not need to repeat the description here except to say that each multiscale 
cross grid has three different mesh sizes, the smallest one being M 3 = 15mm for msgl 
and M 3 = 21mm for msg2. iKrogstad fe Davidson! (|201 lh were careful to design their 



two multiscale cross grids and one conventional grid in such a way that the longitudinal 
integral length-scale of the turbulence at a 2m distance from the grid location is the 
same £q « 23.65mm ± 0.25mm for all three grids. The ratio between £ and the distance 
between the tunnel walls is smaller than 1/75. 

A description of the wind tunnel u sed by IKrogstad fe Davidson! ([201 lh can be found 
in IKrogstad fe Davidson J 2010l 2011 ) . It is much larger and longer than the tunnel used 
by iHurst fe Vassilicosl ( 2007h and the grids were placed in the tunnel contraction, specif- 
ically 1.2m upstream fr om the start of their test sect ion. As a result, the multiscale grid- 
generated turbulence of Krogstad fe Davidson! ( 2011 ) is mo re isotropic and much further 
downstream of the grid than in IHurst fe Vassilicosl (|2007h . Their turbulence measure- 
ments were taken using single and two component hot-wire anemomctry from x ~ 60Z 
till x = 4004 which means 93A/ 3 < x 629M 3 fo r msgl and 67M 3 < x < 446M 3 for 



msg2. This is clearly much further downstream than lHurst fc Vassilicosl ([2007h who could 
not take measurements beyond a distance e qual to 80 times the smallest mesh size of 
their own multiscale cross grids. In the case of Krogstad fc Davidson! ( 2011 ) conventional 
grid (refcred to as eg), 60£o ^ x ^ 400€o corresponds to 40M ^ x ^ 240 M where M is 
the mesh size of the grid. 



The first main conclusion of IKrogstad fc Davidsqn| (1201 lh was that their "r esults are at 
odds" and that their "findings contradict" those of lHurst fc Vassilicosl ([20071 ). This claim 
is factually incorrect not only because the multiscale cross grids in these two studies have 
significant differences (dif ferent blockages, but al s o diffe rences in some other grid-defining 
param eters, see figure 1 in lKrogstad fc Davidson ( 2011 ) and figure 3 in Hurst fc Vassilicoa 
( 20071 ) and related parameters), but also because the regions of the flow were different, 
in fact not even overlapping (if just barely in one case), in terms of multiples of the 
smallest mesh size of the multiscale cross grids. The fact that IHurst fc Vassilicosl ([20071 ) 
measured higher turbulence levels and local Reynolds nu mbers Re\ = uX/v (where A is 



the Taylor microscale and v the kinematic vi scosity) than [K rogsta d fc Davidson! (|201lh 
is quite simply c onsistent with the fact tha t Hurst fc Vassilicos ( 2007h measured much 
closer to the grid. Hurst fc Vassilicosl ( 20071 ) also extracted decay exponents n from their 
data and showed how they vary continuously as the choice of virtual origin xo varies (see 
their figure 10). Choosing xq close to 0, they found exponents n close to 1.2 for their 
multiscale cross grids and close to 1.4 for their conventional grid with large mesh size. 
However, they also proposed a A-based method for choosing xq which returned n « 1.75 
for their multiscale cross grids, n w 2.3 for their conventional grid with high mesh size 
and n « 1.39 for their conventional grid with usual mesh size. These parti c ular A -based 
exponents are indeed higher than those claimed by IKrogstad fc Davidson! (|201lh which 
lie between 1.12 and 1.25 but they were obtained in completely different regions of the 
flow and, most imp ortantly, with differ e nt fitt ing methods which yielded very different 
values of xq . In fact IHurst fc Vassilicosl (|2007l ) did not include these exponents in their 
conclusions because "more extensive checks of large-scale and small-scale isotropy as well 
as homogeneity will be required to fully conclude on the nature of the turbulence decay 
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behind fractal cross grids, in particular in order to assess the viability of our A-based 
method for estimating xo in the cross grid-generated flows" . 

The second main conclusion of iKrogstad fc Davidson is that, in the far-region 

where they measure, their multiscale cross grids and their conventional grid produce 
"virtually identical" turbulence behavior. Furthermore, quoting from their conclusion, 
"Saffman's decay law is reasonably robust, since the energy decay exponents for all three 
grids are close to Saffman's c lassical prediction of n = 6/ 5". In the next two sections 
we use the data published by IKrogstad fe Davidson ( 2011 ) and show that an attentive 
analysis of their data based on Krogstad fe Davidson ( 2010l) leads to very different con- 
clusions. 



2. Decaying homogeneous isotropic turbulence with three different 
initial conditions 



Krogstad fe Davidson! (|201ll ) established that their turbulent flows were reasonably 



homogeneous at x beyond 2m in terms of longitudinal profiles of variances, skewnesses 
and flatnesses of the streamwise fluctuating velocity component. Their centreline mean 
streamwise flow U remains constant to within less than ±0.1% for all three grids from 
x = 2m till about x — 8m, though it deviates a very little bit for msg2 beyond x ~ 
6.5m. As a result, they chose to design their three grids in such a way that they all 
generate turbulence with nearly same longitudinal integr al length-scale £g at x = 2m . 
The positions x = 80£o fall around 1.9m for all the grids. (jKrogstad fe Davidson! (|2011l) 
in fact recorded, and in a few instances used for their analysis, a few measurements at 
closer distances to the grid, i.e. x as small as about 4W -) The longitudinal length-scale 
£ grows as the turbulence moves downstream, but the ratio between I and the distance 
between the tunnel walls remains very small, for example less than about 1/40 at about 
8m from the grid location. 

They also calculated ratios < u 2 > / < u 2 >, < u 2 > / < u 2 z > and u 2 j < u 2 > and 
found small levels of anisotropy "comparable, if not better, than in most other experi- 
ments". In particular, u 2 / < u 2 > hovers between 0.95 and 1.02 throughout the regions 
where they recorded t heir measurements. Hence any anisotropy-rclatcd dependence on 
initial conditions as in Lavoie et all ( 2007 ) can, most probably, be ruled out. 

In figure [Tit we plot < u 2 > jU 2 versus (x — xo)/£o for all three grids as well as fits of 
the data by < u 2 > /U 2 ~ ( x ^" )~ n - The decay exponents n and virtual origins x in 
these fits are estimated simultaneously by direct application of a non-linear least-squares 
regression algorithm ('NLI NFIT' routin e in M ATLAB™). This fitting method is closely 
relate d to the one used by Lavoie et al\ ( 2007 ) and we apply it to nearly the same range 
where IKrogstad fe Davidson ( 201 if ) applied their own fitting methods. Specifically, we 
apply our fit to the range 80£q < x < 330£q which is a range of x fr om about 1.9m to 8m. 
This m eans that, for each grid, we exclude data points obtained bv IKrogstad fe Davidson 
( 201ll ) at values of x smaller than 80£q where according to these authors the turbulence 
is not sufficiently homogeneous, and we also exclude, exactly like IKrogstad fe Davidson 
( 201 lh do, the data points furthest downstream where noise starts to be sign i ficant . 
(Including data points from x » 1.5m (i.e. 60£q) as in IKrogstad fc Davidson ( 2011 ) 
makes little difference as the values of n remain the same to within ±0.01.) 

We give the values of n and xq thus obtained in table [TJ ( method I). These value s agree 
fairly well with the various values of n and x obtained bv IKrogstad fc Davidson ( 2011 ) 
by their three different fitting methods for all three grids except for their value of xo for 
rasgl and their value of n for msgl when they use one of their three fitting methods, the 
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regression method (see their table 1). The values of n which they obtain for msgl with 
their other two fitting methods are close to our value of n for msgl. 

At this point it may be helpful to recall some basic theoretical considerations. Homo- 
geneous turbulence in the wind tunnel decays according to U-^^u 2 = — e where e is the 
turbulent kinetic energy dissipation per unit mass. To obtain and the numerical 

value of n, one needs some more information about u 2 and e. This information usually 
consists of the following three ing redients when the homogeneous turbulence c a n also 
be co nsidered fairly isotropic (see Batchelor fc Townsendl Il948t iBatchelor 19531 : iRottal 
19721 ) : (i) a finite invariant of the von Karman-Howarth equation, (ii) the assumption 
that the decay of large eddies is self-similar and (iii) the empirical assumption that 



A = e£/i 



(2.1) 



remains constant during decay {£ = £(x) is the longitudinal integral length-scale). This 
constancy can be thought of as resulting from the assumed independence of A on turbu- 
le nce intensi t y and Re\. 



Vassilicosl (|201l[ ) proved that there arc four different cases of finite invariants of the 
von Karman-Howarth equation depending on conditions at infinity. A case where no 
known finite invariant exists; a case where the Loitsyansky invariant is the only known 
finite invariant and where self-similar decay of large eddies implies u 2 £ 5 = const during 
decay; a case where only one known finite invariant exists and where self-similar decay 
of large eddies implies u 2 £ m+1 = const with 2 ^ m < 4 (2 ^ ro en s ures t hat the 
spectral tens or does not diverge at zero w avenum ber as stated in IRottal (|l972h . in the 
Appendix of Krogstad fc Davidson! ([201 if ) and in IVassilicod (|201llV )7 and a case where 
two finite invariants exist and where, as a consequence, self-similar decay of large-eddies 
is impossible. 

Using the constancies of A and u 2 £ m+1 , the second and third of these four cases imply 



n = 2(m+ l)/(m + 3) 



(2.2) 



where 2 ^ m ^ 4 and therefore 6/5 ^ n ^ 10/7. There is no known way to rule out the 
first and fourth cases and therefore no known theoretical reason for measured values of 
n to necessarily lie inside the range 6/5 ^ n ^ 10/7. 

Two out of the three present grids have returned values of n which are below 6/5 = 1.2 
(see table Q] under method I). However, this does not imply that the present turbulenc e 
measurements do not fall under the second or third cases identified bv I Vassilicosl (|201ll ). 
Indeed, as lKrogstad fc Davidson! ( 2010l 2011 ) have observed, A varies slowly with x and 
is therefore not strictly constant. If this is so, then (|2.2|) needs to change. 

In figure [2] we plot the values of A obtained by Krogstad fc Davidson! (l201l[ ) for 
thei r three grids as functi o ns of (x — x )/£ where xo is taken from table Q] (method 
I). (jKrogstad fc Davidson! (j 2 1 lh assumed small-scale isotropy and calculated A from 
measurements of < (^f) 2 > using e = 15v < (^f) 2 > and integrations of measured 
longitudinal correlation functions to educe £.) To bring out more clearly the differences 
between grids we in fact plot Aj A\ wher e A-\ is the value of A obtain ed at the smallest 
distance x from each grid. We then follow Krogstad fc Davidson! ( 2010h and fit the power 
law A ~ { x ~ i x ° )~ p in the range 60£o < x < 3304) of this data. These fits are shown in 
figure [5] and the values of p are reported in table [1] 

If A = const is replaced by A ~ (x — xo)~ p then the implication of u 2 £ m+1 = const 
changes from (|2.2I) to 

n= (l-p)2(m + l)/(m + 3) (2.3) 
where 2 ^ m ^ 4. With our estimates of n and p we can now use (|2.3|) to derive values of 
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m for each grid. They are given in table [T] (under method I) and, having now taken into 
account the slight variations of A, they are all between 2 and 4. Similarly, the values of 
n corr = n/(l — p) lie all between 6/5 and 10/7 (see tabic [T] under method I). 

These values of m raise the possi bility that the three decaying nearly homogeneous and 
nearly isotropic turbulent flows of IKrogstad fc Davidson ( 2011 ) may be three different 
instances of the third case identified by IVassilicosI (|201ll ) where only one known finite 



invariant exists and where self-similar decay of large eddies implies u 2 £ m+1 = const 
with 2 ^ m < 4. The Saffman invar iant corresponds to m = 2 but none of the grids 
used by IKrogstad fc Davidson! (|2011l) returns such a value of m. In figure [5Jt we plot 



< u 2 > £ m+1 /(U 2 £ 7 q + ) versus (x — x n )/£ with the values of m given under method 
I in table [TJ for each one of the three different turbulent flows. This figure should be 
compa red with figure [3b which is a reproduction of figure 10 in IKrogstad fc Davidson 



(|201ll ) where they plotted < u 2 x > £ 3 /(U 2 £q) versus {x — xo)/^o, except that we have 



offset the data vertically so as to see more clearly the differences in behavior between each 
grid. Assuming the turbulence i s sufficiently homogeneous a nd isotropic and equally so 
for all three flows (as claimed by Krogstad &: Davidsonll2011 ) , it is clear that the Saffman 
prediction is not satisfied in these flows. Instead, 



2»m+l 



const 



(2.4) 



with m > 2.5 for all grids in the range lOO^o ^ % — %o ^ 400^o- Furthermore, different 
grids give rise to different values of m reaching up to m w 3 with method I (see table [1]). 

In fact there is another way to extract values for n and m from the data (method II), 
and this way gives even better defined invariants and even greater differences between 
the far downstream turbulence decays originating from the conventional grid and the 
multiscalc cross grids. Method II is based on figure lb. This figure is a log-log plot of 
Re\/Re\i versus (x — xq)/Io where Re\\ is the value of Re\ at the smallest distance 
from each grid on this plot and xo is the virtual origin obtained from our nonlinear fit of 
figure la. The first inescapable observation is that the streamwise distributions of Re\ 
are clearly different for the conventional grid and for the multiscalc grids. 

The pow er law form (11.11) implies A 2 ~ (x — xq) in decaying homogeneous isotropic 



turbulence ( Batchelorl (|1953l )). It follows that Re\ ~ {x — xo)*- 1 ")/ 2 , so that a best fit 



of the data in figure lb gives values of n. We apply this power law fit to the very same 
range 801q < x < 330Zo used in method I for our fit of the turbulence intensity data in 
figure la. The values of n thus obtained, the resulting values of m using (|2.3|) and the 
resulting n corr = n/(l + p) arc given in table [TJ under method II. In figure [Jp we use 
these new values of m to plot < u 2 > l m+1 / (U 2 £™ +1 ) versus (x — xq)/£o and find that 
they yield even better defined invariants (|2.4I) than method I (compare with figure [3^.). 
The difference between values of m for conventional grids and values of m for multiscalc 
grids is unmistakable and even greater with method II than with method I. 

We must conclude that the decay of approximately homogeneous turbulence far from 
its initial conditions remains dependent on these initial conditions. The decay exponent 



n and the conserved finite invariant 



both clearly change when the turbulence- 



generating grid is changed. These initial conditions may have to do with the geometry 
of the grids or/and with the inlet Reynolds numbers as the mean speed in the tunnel 
was 13.5m/s when the conventional grid was tested, 14.0m/s when msgl was tested and 
15.5m/s when msg2 was tested. 

In trying to identify the flow-relevant geometrical varia tions from grid to grid, we note 
that the mustiscale grids of IKrogstad fc Davidson! (|201ll ) have three different bar widths 
t\ = 8mm, t% = 4mm and £3 = 2mm as well as three different mesh sizes Mi = 2M2 and 
Mi rts 2M 3 with Mi = 64mm in the case of msgl and Mi = 88mm in the case of msg2 
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30 50 100 300 500 '30 50 100 300 500 



(x-x )/to (x-x a )/l a 

Figure 1. Data and best-fit power laws for: (a) < > /U? nt , an versus (x—xo)/£o (b) Re\/Re\i 
versus (x — xq)/Iq. (•) eg, (A) msgl, (□) msg2. The vertical dashed lines mark the start and 
end of the admissible data range used in the least-squares fits. 



Table 1. Estimation of quantities via least squares fit 

Method I Method II 



Grid p n xo(m) m n corr a n m n corr a 

eg 0.126 1.13 0.23 2.67 1.29 1.90 1.15 2.85 1.32 1.68 

msgl 0.101 1.18 0.28 2.79 1.31 1.14 1.24 3.38 1.37 0.86 

msg2 0.072 1.23 0.33 2.94 1.33 0.62 1.25 3.14 1.35 0.57 



(see their figure 1). Hence the ratio of mesh size to bar width is 11 for msg2 and 8 for 
msgl. It is therefore at least double than M/t = 4 for the conventional grid eg. 

The mesh size determines the distance between the wakes of the bars and the bar 
thickness determines the width of these wakes. Hence the ratio of mesh size to bar 
thickness determines the distance from the grid where the wakes meet and this distance 
increases when we move from eg to msgl and msg2. 

The Reynolds numbers characterising these wakes (calculated as the mean flow speed 
multiplied by the bar thickness and divided by the kinematic viscosity of the air) take 
the values 3.6 x 10 3 in the case of eg; 3 x 10 3 , 1.5 x 10 3 and 7.5 x 10 2 in the case of msgl ; 
and 3.32 x 10 3 , 1.66 x 10 3 and 8.3 x 10 2 in the case of msg2. Unlike conventional grids, 
multiscalc grids impose more than one Reynolds number on the flows they generate and 
a number of different distances from the grid where wakes of different sizes meet. Of 
course the largest wakes are affected by the wakes generated by the smaller ones. But it 
is clear that the turbulence undergoes different generation mechanisms extending over 
different strcamwise distances with different grids. It is indeed remarkable that memory 
of these mechanisms remains in the values of n and m as far downstream as where 



Krogstad fc Davidson! (|201 lh took their measurements. 



3. Different far-field \ow-Re\ turbulent flows 

We obtain ed figures EK.b by taking into account the slow streamwise variation of A as 
suggested by lKrogstad k Davidson] pOlut l. This strcamwise variation can result from the 
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Figure 2. Data and best-fit power laws for A/Ai versus (x — xo)/io- (*) eg, (A) msgl, (d) 
msg2. The vertical dashed lines mark the start and end of the admissible data range used in the 
least-squares fits. 



5 



fc 3 



50 100 150 200 250 300 350 400 450 

(x - x )/e 



50 100 150 200 250 300 350 400 450 

(x - x a )/e 
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(x - X )/i 

Figure 3. Checks of invariant forms via plots of u 2 £ m+1 : (a) m from method I, (b) m from 
method II, (c) m = 2 corresponding to Saffman turbulence. (•) eg, (A) msgl, (D) msg2. For 
improved readability the eg /msgl /msg2 data were vertically offset by (a) [2.0, 1.0, 0.0] x 10~ 4 , 
(b) [1.6, 2.5, 0.0] x 10~ 4 , (c) [2.5, 0.5, -1.5] x 10~ 4 .The left and right vertical dashed lines mark 
the start and end of the data range used to obtain the decay exponents n. This is the range not 
significantly affected by inhomogeneity (to the left) and noise (to the right). 
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Fi gure 4. Re\ dependence of t he normalised energy dissipation rate A (note that A here is 3 A/2 
in IKrogstad fc Davidsonl (|201lf ). for example see their figure 11): (a) eg, (b) msgl, (c) msg2, (d) 
logarithmic axes. (•) eg, (A) msgl, ( n ) msg2. The solid lines are plots of A = const x J?e" 
with a taken from table [TJ method II. 



well-known dependence t hat the dimensionless dissipation ra te A has on Re\ when Re\ 
is below at least 100 (e.g. Burattini. Lavoie fc Antoniall2005[). Indeed, the v alues of Re\ 



characterising the three far- field turbulent flows of IKrogstad fc Davidsonl (|2011 ) range 



between about 90 near x « 604 and 70 at x ps 330£ - Using Re x ~ (x - x ) (1 "™ )/2 and 
A ~ [x — xq)~ p we obtain 

A - Re c l (3.1) 

where 

a = 2p/(n-l). (3.2) 

The values of a implied by this formula on the basis of the exponents p and n obtained in 
the previous section are very different for different grids, ranging from a ~ 1.7 toa« 0.6 
(using n obtained from method II, see table [TJ. Figure |4] confirms how dramatically dif- 
ferent the dependencies of A on Re\ are for the multiscale grids and for the conventional 
grid. Hence, multiscale grids definitely do not "produce almost iden tical results to the 
equivalent classical grids" as claimed by IKrogstad fc Davidson ( 2011 ) . 



Figure |4] also shows that (|3.1|) - (|3.2p give rise to more or less reasonable fits of the data 
thus lending support to the idea that much of the streamwise variation of A comes from 
its dependence on Re\. Increasing values of the dimensionless dissipation rate A with 
increasing Re\ have also been reported in previous works with square bar grids at such 
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rela tively low Reynolds numbers, se e for example fi gure 1 inlBurattini et aLl (12005ft . table 
4 in lComte-Bellot fe Corrsml (|l97ll ) and table 3 in lGad-El-Hak fe CorrsiiTll 19741 ). 



4. Conclusion 

According to the published data in lKrogstad fe Davidson ( 2011 ). multiscale cross grids 
and their equivalent (in terms of Iq) conventional grid can produce very different far- field 
approximately homogeneous isotropic turbulence with wide variations in the dimension- 
less dissipation r ate's dependence on Re \. This would seem to confirm the observation 
already made by iBurattini et all (|2005j ) on the basis of different Re\ dependencies of 
A for different grids, namely that "the geometry of the grid appears to have a per- 
sistent influence in the stre amwise direction up to x/M = 80". In fact the data of 
Krogstad fe Davidson! ( 2011 ) extend this observation to much further distances down- 



stream and to a wider range of grids. 

This data also leads to the conclusi on that the decay of the thre e approximately ho- 
mogeneous isotropic turbulent flows of Krogstad fe Davidson! ( 2011 ) is characterised by 
an invariant quantity u 2 £ m+1 in the region of the flow x ^ 80£o which is the most clearly 
homogeneous. The exponent m is signific antly different from Saffman's m = 2 and ranges 
between 2.7 and 3.4 for the grids used by Krogstad fe Davidson! (120111 ). Their multiscale 



grids return values of m which arc markedly larger than the values of m returned by 
their conventional grid. The streamwise distributions of Re\ and A are also very clearly 
different. 

Fin ally we repeat the remark in our introduction that the data of Krogstad fc Davidson! 
( 201lft do not contradict previous findings on multiscale grids but in fact complement 
them. 



We thank Prof essor Per- Age Krogstad for ki ndly providing us with the post-processed 
data published in iKrogstad fc Davidson! (|201ll ). 
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